MetaboAnalystR 4.0: a unified LC-MS workflow for global metabolomics.

The wide applications of liquid chromatography - mass spectrometry (LC-MS) in untargeted metabolomics demand an easy-to-use, comprehensive computational workflow to support efficient and reproducible data analysis. However, current tools were primarily developed to perform specific tasks in LC-MS based metabolomics data analysis. Here we introduce MetaboAnalystR 4.0 as a streamlined pipeline covering raw spectra processing, compound identification, statistical analysis, and functional interpretation. The key features of MetaboAnalystR 4.0 includes an auto-optimized feature detection and quantification algorithm for LC-MS1 spectra processing, efficient MS2 spectra deconvolution and compound identification for data-dependent or data-independent acquisition, and more accurate functional interpretation through integrated spectral annotation. Comprehensive validation studies using LC-MS1 and MS2 spectra obtained from standards mixtures, dilution series and clinical metabolomics samples have shown its excellent performance across a wide range of common tasks such as peak picking, spectral deconvolution, and compound identification with good computing efficiency. Together with its existing statistical analysis utilities, MetaboAnalystR 4.0 represents a significant step toward a unified, end-to-end workflow for LC-MS based global metabolomics in the open-source R environment.

MetaboAnalyst 6.0: towards a unified platform for metabolomics data processing, analysis and interpretation.

We introduce MetaboAnalyst version 6.0 as a unified platform for processing, analyzing, and interpreting data from targeted as well as untargeted metabolomics studies using liquid chromatography - mass spectrometry (LC-MS). The two main objectives in developing version 6.0 are to support tandem MS (MS2) data processing and annotation, as well as to support the analysis of data from exposomics studies and related experiments. Key features of MetaboAnalyst 6.0 include: (i) a significantly enhanced Spectra Processing module with support for MS2 data and the asari algorithm; (ii) a MS2 Peak Annotation module based on comprehensive MS2 reference databases with fragment-level annotation; (iii) a new Statistical Analysis module dedicated for handling complex study design with multiple factors or phenotypic descriptors; (iv) a Causal Analysis module for estimating metabolite - phenotype causal relations based on two-sample Mendelian randomization, and (v) a Dose-Response Analysis module for benchmark dose calculations. In addition, we have also improved MetaboAnalyst's visualization functions, updated its compound database and metabolite sets, and significantly expanded its pathway analysis support to around 130 species. MetaboAnalyst 6.0 is freely available at https://www.metaboanalyst.ca.

The Native Microbiome Member Chryseobacterium sp. CHNTR56 MYb120 Induces Trehalose Production via a Shift in Central Carbon Metabolism during Early Life in C. elegans.

Aging is the system-wide loss of homeostasis, eventually leading to death. There is growing evidence that the microbiome not only evolves with its aging host, but also directly affects aging via the modulation of metabolites involved in important cellular functions. The widely used model organism C. elegans exhibits high selectivity towards its native microbiome members which confer a range of differential phenotypes and possess varying functional capacities. The ability of one such native microbiome species, Chryseobacterium sp. CHNTR56 MYb120, to improve the lifespan of C. elegans and to promote the production of Vitamin B6 in the co-colonizing member Comamonas sp. 12022 MYb131 are some of its beneficial effects on the worm host. We hypothesize that studying its metabolic influence on the different life stages of the worm could provide further insights into mutualistic interactions. The present work applied LC-MS untargeted metabolomics and isotope labeling to study the impact of the native microbiome member Chryseobacterium sp. CHNTR56 MYb120 on the metabolism of C. elegans. In addition to the upregulation of biosynthesis and detoxification pathway intermediates, we found that Chryseobacterium sp. CHNTR56 MYb120 upregulates the glyoxylate shunt in mid-adult worms which is linked to the upregulation of trehalose, an important metabolite for desiccation tolerance in older worms.

Genomic and phenotypic profiling of Staphylococcus aureus isolates from bovine mastitis for antibiotic resistance and intestinal infectivity.

BACKGROUND: Staphylococcus aureus is one of the prevalent etiological agents of contagious bovine mastitis, causing a significant economic burden on the global dairy industry. Given the emergence of antibiotic resistance (ABR) and possible zoonotic spillovers, S aureus from mastitic cattle pose threat to both veterinary and public health. Therefore, assessment of their ABR status and pathogenic translation in human infection models is crucial. RESULTS: In this study, 43 S. aureus isolates associated with bovine mastitis obtained from four different Canadian provinces (Alberta, Ontario, Quebec, and Atlantic provinces) were tested for ABR and virulence through phenotypic and genotypic profiling. All 43 isolates exhibited crucial virulence characteristics such as hemolysis, and biofilm formation, and six isolates from ST151, ST352, and ST8 categories showed ABR. Genes associated with ABR (tetK, tetM, aac6', norA, norB, lmrS, blaR, blaZ, etc.), toxin production (hla, hlab, lukD, etc.), adherence (fmbA, fnbB, clfA, clfB, icaABCD, etc.), and host immune invasion (spa, sbi, cap, adsA, etc.) were identified by analyzing whole-genome sequences. Although none of the isolates possessed human adaptation genes, both groups of ABR and antibiotic-susceptible isolates demonstrated intracellular invasion, colonization, infection, and death of human intestinal epithelial cells (Caco-2), and Caenorhabditis elegans. Notably, the susceptibilities of S. aureus towards antibiotics such as streptomycin, kanamycin, and ampicillin were altered when the bacteria were internalized in Caco-2 cells and C. elegans. Meanwhile, tetracycline, chloramphenicol, and ceftiofur were comparatively more effective with ≤ 2.5 log10 reductions of intracellular S. aureus. CONCLUSIONS: This study demonstrated the potential of S. aureus isolated from mastitis cows to possess virulence characteristics enabling invasion of intestinal cells thus calling for developing therapeutics capable of targeting drug-resistant intracellular pathogens for effective disease management.

Native Microbiome Members of C. elegans Act Synergistically in Biosynthesis of Pyridoxal 5'-Phosphate.

The roles of the healthy microbiome on the host and the relationships between members of the microbiome remain to be fully characterized. Due to the complexity of the interactions between the mammalian microbiome and its host, the use of model organisms such as the nematode worm Caenorhabditis elegans is a promising strategy to study host-microbiome interactions in vivo, as well as bacterial crosstalk within the host. Previously it was found that native bacterial isolates of the worm, Chryseobacterium sp. CHNTR56 MYb120 and Comamonas sp. 12022 MYb131, possess genomic diversity in the biosynthesis of the active form of vitamin B6, pyridoxal 5'-phosphate (PLP), and contribute to host fitness and lifespan extension. However, the relative contribution of PLP from each isolate, as well as the existence of interbacterial relationships within the worm gut remain to be characterized. In the present work, we investigated the presence and measured the abundance of PLP in the isolates and in the worms grown with the isolates using ultraperformance liquid chromatography tandem-mass spectrometry (UPLC-MS/MS). Our analyses confirmed the presence of PLP in vitro and in vivo. The elevated abundance of PLP in the isolates (which reached statistically significant levels when the two isolates were combined), and within worms grown with the combination of bacterial isolates, compared to control, indicated synergism between the isolates in the production of PLP. Isotope labeling revealed that Comamonas sp. 12022 MYb131 was the main provider of PLP in worms grown with the combination of bacterial isolates. The dominance of this isolate inside the worm was further confirmed by a colonization assay. An untargeted metabolomics analysis of the bacteria showed that the pathways related to cell growth, protein synthesis and lipid synthesis/energy production were regulated in the combination group in comparison with Comamonas sp. 12022 MYb131 alone. Furthermore, glutamine, involved in the de novo synthesis of purine and pyrimidines, was specifically abundant in this group, indicating the potential role of this metabolite in initiating and sustaining bacterial growth. This bacterial crosstalk is suggested to promote the growth of Comamonas sp. 12022 MYb131 in vivo, and synthesis of bacterial metabolites such as PLP in the worm gut.

Ascaris suum Informs Extrasynaptic Volume Transmission in Nematodes.

Neural circuit synaptic connectivities (the connectome) provide the anatomical foundation for our understanding of nematode nervous system function. However, other nonsynaptic routes of communication are known in invertebrates including extrasynaptic volume transmission (EVT), which enables short- and/or long-range communication in the absence of synaptic connections. Although EVT has been highlighted as a facet of Caenorhabditis elegans neurosignaling, no experimental evidence identifies body cavity fluid (pseudocoelomic fluid; PCF) as a vehicle for either neuropeptide or biogenic amine transmission. In the parasitic nematode Ascaris suum, FMRFamide-like peptides encoded on flp-18 potently stimulate female reproductive organs but are expressed in cells that are anatomically distant from the reproductive organ, with no known synaptic connections to this tissue. Here we investigate nonsynaptic neuropeptide signaling in nematodes mediated by the body cavity fluid. Our data show that (i) A. suum PCF (As-PCF) contains a catalog of neuropeptides including FMRFamide-like peptides and neuropeptide-like proteins, (ii) the A. suum FMRFamide-like peptide As-FLP-18A dominates the As-PCF peptidome, (iii) As-PCF potently modulates nematode reproductive muscle function ex vivo, mirroring the effects of synthetic FLP-18 peptides, (iv) As-PCF activates the C. elegans FLP-18 receptors NPR-4 and -5, (v) As-PCF alters C. elegans behavior, and (vi) FLP-18 and FLP-18 receptors display pan-phylum distribution in nematodes. This study provides the first direct experimental evidence to support an extrasynaptic volume route for neuropeptide transmission in nematodes. These data indicate nonsynaptic signaling within the nematode functional connectome and are particularly pertinent to receptor deorphanization approaches underpinning drug discovery programs for nematode pathogens.

The symbiotic relationship between Caenorhabditis elegans and members of its microbiome contributes to worm fitness and lifespan extension.

BACKGROUND: A healthy microbiome influences host physiology through a mutualistic relationship, which can be important for the host to cope with cellular stress by promoting fitness and survival. The mammalian microbiome is highly complex and attributing host phenotypes to a specific member of the microbiome can be difficult. The model organism Caenorhabditis elegans and its native microbiome, discovered recently, can serve as a more tractable, experimental model system to study host-microbiome interactions. In this study, we investigated whether certain members of C. elegans native microbiome would offer a benefit to their host and putative molecular mechanisms using a combination of phenotype screening, omics profiling and functional validation. RESULTS: A total of 16 members of C. elegans microbiome were screened under chemically-induced toxicity. Worms grown with Chryseobacterium sp. CHNTR56 MYb120 or Comamonas sp. 12022 MYb131, were most resistant to oxidative chemical stress (SiO2 nanoparticles and juglone), as measured by progeny output. Further investigation showed that Chryseobacterium sp. CHNTR56 positively influenced the worm's lifespan, whereas the combination of both isolates had a synergistic effect. RNAseq analysis of young adult worms, grown with either isolate, revealed the enrichment of cellular detoxification mechanisms (glutathione metabolism, drug metabolism and metabolism of xenobiotics) and signaling pathways (TGF-beta and Wnt signaling pathways). Upregulation of cysteine synthases (cysl genes) in the worms, associated with glutathione metabolism, was also observed. Nanopore sequencing uncovered that the genomes of the two isolates have evolved to favor the specific route of the de novo synthesis pathway of vitamin B6 (cofactor of cysl enzymes) through serC or pdxA2 homologs. Finally, co-culture with vitamin B6 extended worm lifespan. CONCLUSIONS: In summary, our study indicates that certain colonizing members of C. elegans have genomic diversity in vitamin B6 synthesis and promote host fitness and lifespan extension. The regulation of host cellular detoxification genes (i.e. gst) along with cysl genes at the transcriptome level and the bacterium-specific vitamin B6 synthesis mechanism at the genome level are in an agreement with enhanced host glutathione-based cellular detoxification due to this interspecies relationship. C. elegans is therefore a promising alternative model to study host-microbiome interactions in host fitness and lifespan.

Comprehensive phenotyping and transcriptome profiling to study nanotoxicity in C. elegans.

Engineered nanoparticles are used at an increasing rate in both industry and medicine without fully understanding their impact on health and environment. The nematode Caenorhabditis elegans is a suitable model to study the toxic effects of nanoparticles as it is amenable to comprehensive phenotyping, such as locomotion, growth, neurotoxicity and reproduction. In this study, we systematically evaluated the effects of silver (Ag) and five metal oxide nanoparticles: SiO2, CeO2, CuO, Al2O3 and TiO2. The results showed that Ag and SiO2 exposures had the most toxic effects on locomotion velocity, growth and reproduction, whereas CeO2, Al2O3 and CuO exposures were mostly neurotoxic. We further performed RNAseq to compare the gene expression profiles underlying Ag and SiO2toxicities. Gene set enrichment analyses revealed that exposures to Ag and SiO2consistently downregulated several biological processes (regulations in locomotion, reproductive process and cell growth) and pathways (neuroactive ligand-receptor interaction, wnt and MAPK signaling, etc.), with opposite effects on genes involved in innate immunity. Our results contribute to mechanistic insights into toxicity of Ag and SiO2 nanoparticles and demonstrated that C. elegans as a valuable model for nanotoxicity assessment.

T1000: a reduced gene set prioritized for toxicogenomic studies.

There is growing interest within regulatory agencies and toxicological research communities to develop, test, and apply new approaches, such as toxicogenomics, to more efficiently evaluate chemical hazards. Given the complexity of analyzing thousands of genes simultaneously, there is a need to identify reduced gene sets. Though several gene sets have been defined for toxicological applications, few of these were purposefully derived using toxicogenomics data. Here, we developed and applied a systematic approach to identify 1,000 genes (called Toxicogenomics-1000 or T1000) highly responsive to chemical exposures. First, a co-expression network of 11,210 genes was built by leveraging microarray data from the Open TG-GATEs program. This network was then re-weighted based on prior knowledge of their biological (KEGG, MSigDB) and toxicological (CTD) relevance. Finally, weighted correlation network analysis was applied to identify 258 gene clusters. T1000 was defined by selecting genes from each cluster that were most associated with outcome measures. For model evaluation, we compared the performance of T1000 to that of other gene sets (L1000, S1500, Genes selected by Limma, and random set) using two external datasets based on the rat model. Additionally, a smaller (T384) and a larger version (T1500) of T1000 were used for dose-response modeling to test the effect of gene set size. Our findings demonstrated that the T1000 gene set is predictive of apical outcomes across a range of conditions (e.g., in vitro and in vivo, dose-response, multiple species, tissues, and chemicals), and generally performs as well, or better than other gene sets available.

PmrC (EptA) and CptA Negatively Affect Outer Membrane Vesicle Production in Citrobacter rodentium.

Outer membrane vesicles (OMVs) are naturally produced by Gram-negative bacteria by a bulging of the outer membrane (OM) and subsequent release into the environment. By serving as vehicles for various cargos, including proteins, nucleic acids and small metabolites, OMVs are central to interbacterial interactions and both symbiotic and pathogenic host bacterial interactions. However, despite their importance, the mechanism of OMV formation remains unclear. Recent evidence indicates that covalent modifications of lipopolysaccharides (LPS) influence OMV biogenesis. Several enteric bacteria modify LPS with phosphoethanolamine (pEtN) using the iron-regulated PmrC (EptA) and CptA pEtN transferases. In wild-type Citrobacter rodentium, the presence of increasing subtoxic concentrations of iron was found to stimulate OMV production 4- to 9-fold above baseline. C. rodentium uses the two-component system PmrAB to sense and adapt to environmental iron. Compared to the wild type, the C. rodentium ΔpmrAB strain exhibited heightened OMV production at similar iron concentrations. PmrAB regulates transcription of pmrC (also known as eptA) and cptA OMV production in strains lacking either pmrC (eptA) or cptA was similarly increased in comparison to that of the wild type. Importantly, plasmid complementation of C. rodentium strains with either pmrC (eptA) or cptA resulted in a drastic inhibition of OMV production. Finally, we showed that ß-lactamase and CroP, two enzymes found in the C. rodentium periplasm and outer membrane (OM), respectively, are associated with OMVs. These data suggest a novel mechanism by which C. rodentium and possibly other Gram-negative bacteria can negatively affect OMV production through the PmrAB-regulated genes pmrC (eptA) and cptAIMPORTANCE Although OMVs secreted by Gram-negative bacteria fulfill multiple functions, the molecular mechanism of OMV biogenesis remains ill defined. Our group has previously shown that PmrC (also known as EptA) and CptA maintain OM integrity and provide resistance to iron toxicity and antibiotics in the murine pathogen Citrobacter rodentium In several enteric bacteria, these proteins modify the lipid A and core regions of lipopolysaccharide with phosphoethanolamine moieties. Here, we show that these proteins also repress OMV production in response to environmental iron in C. rodentium These data support the emerging understanding that lipopolysaccharide modifications are important regulators of OMV biogenesis in Gram-negative bacteria.

Absence of PmrAB-mediated phosphoethanolamine modifications of Citrobacter rodentium lipopolysaccharide affects outer membrane integrity.

The PmrAB two-component system of enterobacteria regulates a number of genes whose protein products modify lipopolysaccharide (LPS). The LPS is modified during transport to the bacterial outer membrane (OM). A subset of PmrAB-mediated LPS modifications consists of the addition of phosphoethanolamine (pEtN) to lipid A by PmrC and to the core by CptA. In Salmonella enterica, pEtN modifications have been associated with resistance to polymyxin B and to excess iron. To investigate putative functions of pEtN modifications in Citrobacter rodentium, ΔpmrAB, ΔpmrC, ΔcptA, and ΔpmrC ΔcptA deletion mutants were constructed. Compared to the wild type, most mutant strains were found to be more susceptible to antibiotics that must diffuse across the LPS layer of the OM. All mutant strains also showed increased influx rates of ethidium dye across their OM, suggesting that PmrAB-regulated pEtN modifications affect OM permeability. This was confirmed by increased partitioning of the fluorescent dye 1-N-phenylnaphthylamine (NPN) into the OM phospholipid layer of the mutant strains. In addition, substantial release of periplasmic ß-lactamase was observed for the ΔpmrAB and ΔpmrC ΔcptA strains, indicating a loss of OM integrity. This study attributes a new role for PmrAB-mediated pEtN LPS modifications in the maintenance of C. rodentium OM integrity.

Activity of novel nicotinic anthelmintics in cut preparations of Caenorhabditis elegans.

The free-living nematode Caenorhabditis elegans is a useful model for studying the pharmacology of anthelmintics. Currently approved anthelmintics have various mechanisms of action, including activity at nematode nicotinic acetylcholine receptors (nAChRs). Classical anthelmintic agonists of these receptors (nicotine, levamisole, pyrantel and bephenium) caused intact specimens of C. elegans to undergo contracted paralysis. The nAChR antagonist mecamylamine paralysed intact worms and blocked the actions of the agonists. The time to onset of effects of these drugs was enhanced when worms bisected between the mid- and anterior-portions were tested. The novel anthelmintic nAChR antagonist derquantel (2-desoxoparaherquamide, 2-DOPH) was weakly active in intact specimens of C. elegans at concentrations of 50 µM over several days. No antagonism of the nAChR agonists was observed with this drug in intact worms. However, derquantel had direct and marked effects on motility in cut worms and blocked the effects of nAChR agonists in this preparation. A representative of the new amino-acetonitrile derivative (AAD) class of nAChR agonists was not antagonised by derquantel in cut C. elegans, suggesting that these two anthelmintics may not demonstrate unfavourable drug-drug interactions at the receptor level if used to treat livestock infected with parasitic nematodes. The permeability properties of the C. elegans cuticle may be more restrictive than those of adult parasites, calling into question primary anthelmintic screening strategies that rely on this model organism.

An outer membrane protease of the omptin family prevents activation of the Citrobacter rodentium PhoPQ two-component system by antimicrobial peptides.

The PhoPQ two-component system of the intracellular pathogen Salmonella enterica senses and controls resistance to alpha-helical antimicrobial peptides (AMPs) by regulating covalent modifications of lipid A. A homologue of the phoPQ operon was found in the genome of the murine enteric extracellular pathogen, Citrobacter rodentium. Here we report that C. rodentium PhoPQ was apparently unable to mediate activation of target genes in the presence of alpha-helical AMPs. However, these AMPs activated C. rodentium PhoPQ expressed in a S. entericaDeltaphoPQ mutant. Analysis of the outer membrane (OM) fractions of the C. rodentium wild-type and DeltaphoPQ strains led to the identification of an omptin family protease (CroP) that was absent in DeltaphoPQ. Deletion of croP in C. rodentium resulted in higher susceptibility to alpha-helical AMPs, indicating a direct role of CroP in AMP resistance. CroP greatly contributed to the protection of the OM from AMP damage by actively degrading alpha-helical AMPs before they reach the periplasmic space. Accordingly, transcriptional activation of PhoP-regulated genes by alpha-helical AMPs was restored in the DeltacroP mutant. This study shows that resistance to alpha-helical AMPs by the extracellular pathogen C. rodentium relies primarily on the CroP OM protease.

The prpZ gene cluster encoding eukaryotic-type Ser/Thr protein kinases and phosphatases is repressed by oxidative stress and involved in Salmonella enterica serovar Typhi survival in human macrophages.

The prpZ gene cluster consists of three ORFs coding for proteins with homology to eukaryotic-type Ser/Thr protein phosphatases 2C (prpZ) and Ser/Thr protein kinases (prkY and prkX). This cluster is present in the sequenced genomes of Salmonella enterica serovar Typhi (S. Typhi) strains Ty2 and CT18. This study investigated the genetic organization of this gene cluster, its regulation and its putative involvement in virulence. The three genes are transcribed as a polycistronic mRNA as demonstrated by reverse transcriptase (RT)-PCR. Analysis of a prpZ::lacZ transcriptional fusion showed that the prpZ locus is expressed throughout the growth phase. LacZ activity and real-time RT-PCR showed that transcription of the mRNA is negatively regulated upon exposure of cells to HOCl and, to a lesser extent, hydrogen peroxide. A deletion mutant of the prpZ gene cluster showed a significantly lower level of survival than the parental strain Ty2 in human macrophages at 48 h postinfection. Together these data suggest that prpZ, prkY and prkX are virulence genes that may be part of a signaling pathway controlling long-term survival of S. Typhi in host cells.